Genome organization across scales: mechanistic insights from in vitro reconstitution studies.

Eukaryotic genomes are compacted and organized into distinct three-dimensional (3D) structures, which range from small-scale nucleosome arrays to large-scale chromatin domains. These chromatin structures play an important role in the regulation of transcription and other nuclear processes. The molecular mechanisms that drive the formation of chromatin structures across scales and the relationship between chromatin structure and function remain incompletely understood. Because the processes involved are complex and interconnected, it is often challenging to dissect the underlying principles in the nuclear environment. Therefore, in vitro reconstitution systems provide a valuable approach to gain insight into the molecular mechanisms by which chromatin structures are formed and to determine the cause-consequence relationships between the processes involved. In this review, we give an overview of in vitro approaches that have been used to study chromatin structures across scales and how they have increased our understanding of the formation and function of these structures. We start by discussing in vitro studies that have given insight into the mechanisms of nucleosome positioning. Next, we discuss recent efforts to reconstitute larger-scale chromatin domains and loops and the resulting insights into the principles of genome organization. We conclude with an outlook on potential future applications of chromatin reconstitution systems and how they may contribute to answering open questions concerning chromatin architecture.

In vitro reconstitution of chromatin domains shows a role for nucleosome positioning in 3D genome organization.

Eukaryotic genomes are organized into chromatin domains. The molecular mechanisms driving the formation of these domains are difficult to dissect in vivo and remain poorly understood. Here we reconstitute Saccharomyces cerevisiae chromatin in vitro and determine its 3D organization at subnucleosome resolution by micrococcal nuclease-based chromosome conformation capture and molecular dynamics simulations. We show that regularly spaced and phased nucleosome arrays form chromatin domains in vitro that resemble domains in vivo. This demonstrates that neither loop extrusion nor transcription is required for basic domain formation in yeast. In addition, we find that the boundaries of reconstituted domains correspond to nucleosome-free regions and that insulation strength scales with their width. Finally, we show that domain compaction depends on nucleosome linker length, with longer linkers forming more compact structures. Together, our results demonstrate that regular nucleosome positioning is important for the formation of chromatin domains and provide a proof-of-principle for bottom-up 3D genome studies.

The Mediator complex regulates enhancer-promoter interactions.

Enhancer-mediated gene activation generally requires physical proximity between enhancers and their target gene promoters. However, the molecular mechanisms by which interactions between enhancers and promoters are formed are not well understood. Here, we investigate the function of the Mediator complex in the regulation of enhancer-promoter interactions, by combining rapid protein depletion and high-resolution MNase-based chromosome conformation capture approaches. We show that depletion of Mediator leads to reduced enhancer-promoter interaction frequencies, which are associated with a strong decrease in gene expression. In addition, we find increased interactions between CTCF-binding sites upon Mediator depletion. These changes in chromatin architecture are associated with a redistribution of the Cohesin complex on chromatin and a reduction in Cohesin occupancy at enhancers. Together, our results indicate that the Mediator and Cohesin complexes contribute to enhancer-promoter interactions and provide insights into the molecular mechanisms by which communication between enhancers and promoters is regulated.

Assessment of Multiway Interactions with Tri-C.

Tri-C is a chromosome conformation capture (3C) approach that can efficiently identify multiway chromatin interactions with viewpoints of interest. As opposed to pair-wise interactions identified in methods such as Hi-C, 4C, and Capture-C, the detection of multiway interactions allows researchers to investigate how multiple cis-regulatory elements interact together in higher-order structures in single nuclei and address questions regarding structural cooperation between these elements. Here, we describe the procedure for designing and performing a Tri-C experiment.

Analysis of sub-kilobase chromatin topology reveals nano-scale regulatory interactions with variable dependence on cohesin and CTCF.

Enhancers and promoters predominantly interact within large-scale topologically associating domains (TADs), which are formed by loop extrusion mediated by cohesin and CTCF. However, it is unclear whether complex chromatin structures exist at sub-kilobase-scale and to what extent fine-scale regulatory interactions depend on loop extrusion. To address these questions, we present an MNase-based chromosome conformation capture (3C) approach, which has enabled us to generate the most detailed local interaction data to date (20 bp resolution) and precisely investigate the effects of cohesin and CTCF depletion on chromatin architecture. Our data reveal that cis-regulatory elements have distinct internal nano-scale structures, within which local insulation is dependent on CTCF, but which are independent of cohesin. In contrast, we find that depletion of cohesin causes a subtle reduction in longer-range enhancer-promoter interactions and that CTCF depletion can cause rewiring of regulatory contacts. Together, our data show that loop extrusion is not essential for enhancer-promoter interactions, but contributes to their robustness and specificity and to precise regulation of gene expression.

Capture-C: a modular and flexible approach for high-resolution chromosome conformation capture.

Chromosome conformation capture (3C) methods measure the spatial proximity between DNA elements in the cell nucleus. Many methods have been developed to sample 3C material, including the Capture-C family of protocols. Capture-C methods use oligonucleotides to enrich for interactions of interest from sequencing-ready 3C libraries. This approach is modular and has been adapted and optimized to work for sampling of disperse DNA elements (NuTi Capture-C), including from low cell inputs (LI Capture-C), as well as to generate Hi-C like maps for specific regions of interest (Tiled-C) and to interrogate multiway interactions (Tri-C). We present the design, experimental protocol and analysis pipeline for NuTi Capture-C in addition to the variations for generation of LI Capture-C, Tiled-C and Tri-C data. The entire procedure can be performed in 3 weeks and requires standard molecular biology skills and equipment, access to a next-generation sequencing platform, and basic bioinformatic skills. Implemented with other sequencing technologies, these methods can be used to identify regulatory interactions and to compare the structural organization of the genome in different cell types and genetic models.

Dynamic Runx1 chromatin boundaries affect gene expression in hematopoietic development.

The transcription factor RUNX1 is a critical regulator of developmental hematopoiesis and is frequently disrupted in leukemia. Runx1 is a large, complex gene that is expressed from two alternative promoters under the spatiotemporal control of multiple hematopoietic enhancers. To dissect the dynamic regulation of Runx1 in hematopoietic development, we analyzed its three-dimensional chromatin conformation in mouse embryonic stem cell (ESC) differentiation cultures. Runx1 resides in a 1.1 Mb topologically associating domain (TAD) demarcated by convergent CTCF motifs. As ESCs differentiate to mesoderm, chromatin accessibility, Runx1 enhancer-promoter (E-P) interactions, and CTCF-CTCF interactions increase in the TAD, along with initiation of Runx1 expression from the P2 promoter. Differentiation to hematopoietic progenitor cells is associated with the formation of tissue-specific sub-TADs over Runx1, a shift in E-P interactions, P1 promoter demethylation, and robust expression from both Runx1 promoters. Deletion of promoter-proximal CTCF sites at the sub-TAD boundaries has no obvious effects on E-P interactions but leads to partial loss of domain structure, mildly affects gene expression, and delays hematopoietic development. Together, our analysis of gene regulation at a large multi-promoter developmental gene reveals that dynamic sub-TAD chromatin boundaries play a role in establishing TAD structure and coordinated gene expression.

Defining genome architecture at base-pair resolution.

In higher eukaryotes, many genes are regulated by enhancers that are 104-106 base pairs (bp) away from the promoter. Enhancers contain transcription-factor-binding sites (which are typically around 7-22 bp), and physical contact between the promoters and enhancers is thought to be required to modulate gene expression. Although chromatin architecture has been mapped extensively at resolutions of 1 kilobase and above; it has not been possible to define physical contacts at the scale of the proteins that determine gene expression. Here we define these interactions in detail using a chromosome conformation capture method (Micro-Capture-C) that enables the physical contacts between different classes of regulatory elements to be determined at base-pair resolution. We find that highly punctate contacts occur between enhancers, promoters and CCCTC-binding factor (CTCF) sites and we show that transcription factors have an important role in the maintenance of the contacts between enhancers and promoters. Our data show that interactions between CTCF sites are increased when active promoters and enhancers are located within the intervening chromatin. This supports a model in which chromatin loop extrusion1 is dependent on cohesin loading at active promoters and enhancers, which explains the formation of tissue-specific chromatin domains without changes in CTCF binding.

Enhancers predominantly regulate gene expression during differentiation via transcription initiation.

Gene transcription occurs via a cycle of linked events, including initiation, promoter-proximal pausing, and elongation of RNA polymerase II (Pol II). A key question is how transcriptional enhancers influence these events to control gene expression. Here, we present an approach that evaluates the level and change in promoter-proximal transcription (initiation and pausing) in the context of differential gene expression, genome-wide. This combinatorial approach shows that in primary cells, control of gene expression during differentiation is achieved predominantly via changes in transcription initiation rather than via release of Pol II pausing. Using genetically engineered mouse models, deleted for functionally validated enhancers of the α- and ß-globin loci, we confirm that these elements regulate Pol II recruitment and/or initiation to modulate gene expression. Together, our data show that gene expression during differentiation is regulated predominantly at the level of initiation and that enhancers are key effectors of this process.

High-resolution targeted 3C interrogation of cis-regulatory element organization at genome-wide scale.

Chromosome conformation capture (3C) provides an adaptable tool for studying diverse biological questions. Current 3C methods generally provide either low-resolution interaction profiles across the entire genome, or high-resolution interaction profiles at limited numbers of loci. Due to technical limitations, generation of reproducible high-resolution interaction profiles has not been achieved at genome-wide scale. Here, to overcome this barrier, we systematically test each step of 3C and report two improvements over current methods. We show that up to 30% of reporter events generated using the popular in situ 3C method arise from ligations between two individual nuclei, but this noise can be almost entirely eliminated by isolating intact nuclei after ligation. Using Nuclear-Titrated Capture-C, we generate reproducible high-resolution genome-wide 3C interaction profiles by targeting 8055 gene promoters in erythroid cells. By pairing high-resolution 3C interaction calls with nascent gene expression we interrogate the role of promoter hubs and super-enhancers in gene regulation.

The relationship between genome structure and function.

Precise patterns of gene expression in metazoans are controlled by three classes of regulatory elements: promoters, enhancers and boundary elements. During differentiation and development, these elements form specific interactions in dynamic higher-order chromatin structures. However, the relationship between genome structure and its function in gene regulation is not completely understood. Here we review recent progress in this field and discuss whether genome structure plays an instructive role in regulating gene expression or is a reflection of the activity of the regulatory elements of the genome.

The mouse alpha-globin cluster: a paradigm for studying genome regulation and organization.

The mammalian globin gene clusters provide a paradigm for studying the relationship between genome structure and function. As blood stem cells undergo lineage specification and differentiation to form red blood cells, the chromatin structure and expression of the α-globin cluster change. The gradual activation of the α-globin genes in well-defined cell populations has enabled investigation of the structural and functional roles of its enhancers, promoters and boundary elements. Recent studies of gene regulatory processes involving these elements at the mouse α-globin cluster have brought new insights into the general principles underlying the three-dimensional structure of the genome and its relationship to gene expression throughout time.

DeepC: predicting 3D genome folding using megabase-scale transfer learning.

Predicting the impact of noncoding genetic variation requires interpreting it in the context of three-dimensional genome architecture. We have developed deepC, a transfer-learning-based deep neural network that accurately predicts genome folding from megabase-scale DNA sequence. DeepC predicts domain boundaries at high resolution, learns the sequence determinants of genome folding and predicts the impact of both large-scale structural and single base-pair variations.

Dynamics of the 4D genome during in vivo lineage specification and differentiation.

Mammalian gene expression patterns are controlled by regulatory elements, which interact within topologically associating domains (TADs). The relationship between activation of regulatory elements, formation of structural chromatin interactions and gene expression during development is unclear. Here, we present Tiled-C, a low-input chromosome conformation capture (3C) technique. We use this approach to study chromatin architecture at high spatial and temporal resolution through in vivo mouse erythroid differentiation. Integrated analysis of chromatin accessibility and single-cell expression data shows that regulatory elements gradually become accessible within pre-existing TADs during early differentiation. This is followed by structural re-organization within the TAD and formation of specific contacts between enhancers and promoters. Our high-resolution data show that these enhancer-promoter interactions are not established prior to gene expression, but formed gradually during differentiation, concomitant with progressive upregulation of gene activity. Together, these results provide new insight into the close, interdependent relationship between chromatin architecture and gene regulation during development.

A Dynamic Folded Hairpin Conformation Is Associated with α-Globin Activation in Erythroid Cells.

We investigate the three-dimensional (3D) conformations of the α-globin locus at the single-allele level in murine embryonic stem cells (ESCs) and erythroid cells, combining polymer physics models and high-resolution Capture-C data. Model predictions are validated against independent fluorescence in situ hybridization (FISH) data measuring pairwise distances, and Tri-C data identifying three-way contacts. The architecture is rearranged during the transition from ESCs to erythroid cells, associated with the activation of the globin genes. We find that in ESCs, the spatial organization conforms to a highly intermingled 3D structure involving non-specific contacts, whereas in erythroid cells the α-globin genes and their enhancers form a self-contained domain, arranged in a folded hairpin conformation, separated from intermingling flanking regions by a thermodynamic mechanism of micro-phase separation. The flanking regions are rich in convergent CTCF sites, which only marginally participate in the erythroid-specific gene-enhancer contacts, suggesting that beyond the interaction of CTCF sites, multiple molecular mechanisms cooperate to form an interacting domain.

A revised model for promoter competition based on multi-way chromatin interactions at the α-globin locus.

Specific communication between gene promoters and enhancers is critical for accurate regulation of gene expression. However, it remains unclear how specific interactions between multiple regulatory elements contained within a single chromatin domain are coordinated. Recent technological advances which can detect multi-way chromatin interactions at single alleles can provide insights into how multiple regulatory elements cooperate or compete for transcriptional activation. Here, we use such an approach to investigate how interactions of the α-globin enhancers are distributed between multiple promoters in a mouse model in which the α-globin domain is extended to include several additional genes. Our data show that gene promoters do not form mutually exclusive interactions with enhancers, but all interact simultaneously in a single complex. These findings suggest that promoters do not structurally compete for interactions with enhancers, but form a regulatory hub structure, which is consistent with recent models of transcriptional activation occurring in non-membrane bound nuclear compartments.

DOT1L inhibition reveals a distinct subset of enhancers dependent on H3K79 methylation.

Enhancer elements are a key regulatory feature of many important genes. Several general features including the presence of specific histone modifications are used to demarcate potentially active enhancers. Here we reveal that putative enhancers marked with H3 lysine 79 (H3K79) di or trimethylation (me2/3) (which we name H3K79me2/3 enhancer elements or KEEs) can be found in multiple cell types. Mixed lineage leukemia gene (MLL) rearrangements (MLL-r) such as MLL-AF4 are a major cause of incurable acute lymphoblastic leukemias (ALL). Using the DOT1L inhibitor EPZ-5676 in MLL-AF4 leukemia cells, we show that H3K79me2/3 is required for maintaining chromatin accessibility, histone acetylation and transcription factor binding specifically at KEEs but not non-KEE enhancers. We go on to show that H3K79me2/3 is essential for maintaining enhancer-promoter interactions at a subset of KEEs. Together, these data implicate H3K79me2/3 as having a functional role at a subset of active enhancers in MLL-AF4 leukemia cells.